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Characterization and phylogenetic analysis of е mitochondrial 128 rRNA 


gene sequences from nine species of silk-producing insects 
LIU Yan-Qun'?* 1] JIN Xiang-Dong[] QIN Li'[] LU Cheng[] XIANG Zhong-Huaif] 1. College of Bioscience 
and Biotechnology[] Shenyang Agricultural University[] Shenyang 110161[] China[] 2. The Key Sericultural 
Laboratory of Agricultural Ministry[] Southwest University[] Chongqing 400716[] China[] 3. Sericultural Institute 
of Jilin Province[] Yongji[] Jilin 132200[] China[] 
Abstract[] The nucleotide composition and substitution pattern of the partial mitochondrial 128 rRNA gene from 
nine species of silk-producing insects were investigated. The phylogenetic relationships among the nine species 
were then assessed by the methods of unweighted pair group method using arithmetic averaged] UPGMA[TI 
neighbor joining] NJOO minimum evolution] МЕП and maximum parsimony[] МРП. The 128 rRNA gene 
sequences of 427 bp of the wild Antheraea pernyi are identical to the semi-domestic A. perny[] Yuzao No. 
10. The genera Antheraeal] Samial] and Bombyx show obviously different nucleotide composition[] i[] the AT 
content of Bombycidae family is higher than that of Saturniidae[] 10 the nucleotide of Antheraea and Samia 
П Satumiidae[] contains more Ts Шап As[] while more As than Ts exist in Bombyx[] Bombycidae[]. By analyzing 
the nucleotide substitution pattern[] more transitions than transvertions were found among these species except 
Bombyx . All of the phylogenetic trees showed monophyly for the three genera[] and the confidence limits were 
measured more than 96% for the clades. However[] the phylogenetic tree based on UPGMA method supported 
the ancestral position of A. assama within Antheraea species[] whereas all the trees using NJL] МЕП and MP 


methods backed up the ancestral position of A. mylitta[] suggesting that the phylogenetic relationship among 
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Antheraea members need to be further studied with more molecular markers . 


Key words[] Silk-producing insects Bombycidae[] Saturniidae] 12S rRNA | gene[] nucleotide composition[] 


substitution pattern[] phylogenetic relationship 
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Table 1 Тһе nine species of the silk-producing insects used in this study 
ПП ПП 00000 0 00000 
Species Code Variety or source Silk GenBank accession no. 
00 0 0 Satumiidae 
ах a EET En China ПП Oak AY242296 
e Wn а : Eun аа China ПП ow EU 
Ap3 ПО O Unknown ПО Oak AY037834 
OOOD Antheraea yamamai Ayl [] [] Unknown П 0 Tensan AY037835 
Ay2 OO Japan П 0 Тепзап AF029064 
UD UL Antheraea mylitta Am ПП India ПО D Tasar AY037833 
UD DU Antheraea roylei Ar ПО India ПО Oak AY037832 
000D Antheraea proylei Apr ПО India ПО Oak AY037831 
ПО 0 Antheraea assama Aa ПО India 0O00 Muga AY037830 
OOO Samia cynthia ricini Scr ПО India ПОО Ei AY037829 
П D 0 Bombycidae 
O 0 Bombyx mori Bm1 Xiafang[T] [] China О O Mulberry AY048187 
Bm2 Aojuku[T] П Japan О O Mulberry AB083339 
Bm3 Backokjam[T] [] Korea u Mulberry AF149768 
Bm4 C-108[T] O China О O Mulberry AB070264 
Bm5 П [] Unknown О O Mulberry AF029065 
О 0 0 Bombyx mandarina Bmanl U D [ПП Shangdongf] China - DQ263226 
Bman2 О. 0 Japan - AB070263 
Bman3 П O Unknown - AF029066 
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Table 2 Nucleotide composition and skewness of the 125 rRNA gene sequences for silk-producing insects 








LengtH] bp A% G% T% C% A+T % AT skew GC skew 
UD [D Antheraea 
Apl/2 427 38.4 6.6 44.0 11.0 82.4 — 0.068 — 0.250 
Ap3 428 38.6 6.5 43.9 11.0 82.5 - 0.064 — 0.257 
Ayl 431 39.0 6.7 43.2 11.1 82.2 -0.051 — 0.247 
Ay2 430 38.8 6.7 44.0 10.5 82.8 — 0.063 = 0.221 
Am 429 38.5 6.5 44.5 10.5 83.0 = 0.072 — 0.235 
Ar 419 37.9 7.2 44.2 10.7 82.1 = 0.077 -0.196 
Apr 425 39.5 6.4 43.3 10.8 82.8 — 0.046 = 0.256 
Aa 432 38.4 6.7 44.4 10.4 82.8 – 0.072 -0.216 
O 0 Average 38.6 6.7 43.9 10.8 82.5 = 0.064 – 0.234 
ПО D Samia 
Scr 435 39.8 6.9 42.1 11.3 81.9 – 0.028 — 0.243 
UD 0 Bombyx 
Bml 431 42.0 6.3 42.0 9.7 84.0 0 — 0.213 
Bm2 431 42.2 6.0 42.0 9.7 84.2 0.002 – 0.234 
Bm3 431 42.2 6.0 42.0 9.7 84.2 0.002 – 0.234 
Bm4 431 42.2 6.0 42.0 9.7 84.2 0.002 – 0.234 
Bm5 431 42.5 6.0 41.8 9.7 84.3 0.008 — 0.236 
Bmanl 432 42.1 6.0 42.1 9.7 84.2 0 — 0.234 
Bman2 431 42.5 6.0 42.0 9.5 84.5 0.006 = 0.226 
Bman3 432 42.4 6.0 41.9 9.7 84.3 0.006 — 0.236 
O 0 Average 42.3 6.0 42.0 9.7 84.3 0.004 — 0.236 
ПП] Dm 430 37.7 9.8 38.8 13.7 76.5 -0.014 -0.166 


DM[I[] O Drosophila melanogaster. П ЗП The same for Table 3. 


11111111 1111111122 2222222222 222222233 33333333333 3444444444 
222222233 3334606777 8801122233 3334444801 1111111222 344445503 34577777888 9000001122 
5023456801 2366379018 2653446701 3480134480 1345689347 815688995 98701789349 7234580123 





















































Ар1/2 TT-AATTAAT ATAGTTATTT TTCCTAATCA TTTAC-AATT ССТТТТАТТА ТСАААТТТОС ТТСА-—ТТА TATATA-T-T 
Арз A ose eet ОЛЕ ОЛОК eti ode CE rne 
Арт AAT sis Аб I SIS es TT TEE OH oue supp ed ids АЙДАДЫК an Бу MIS STR ini uds 
Ár ААБ: „убы ек АКУ НЫКЫ e ex ТТЛ. Tessa Оаа Ex LT 
Ау1 АТ о spiri paare Buren НЕ, TG. saisis es адир t pe SUR Á A.G 
Ay2 Aon зреет Ар ТАТ TM А к н Бо e Terawis LU 
Am Е T so eae ees ps Жылы АТТАТ.А. Т......... ы ——... AT... T.T- 
Аа АРАТА Т, лене RS wea Салаа ТТАТ, А, TTA...... ГИТЕ Секси» асы К. i uu 
Scr Kui Tus [o A.T С.АТТ-. А, ТТ..А,ТА,. АТ....А,.. AG., АТТС, С АСАТСТА. ТА 
Вт1 КИ Aust АТАА. ААТА. ТТ. А. AACT. А. САС ТТ. AAATAA, ATTT... AAA . NTTTAA, А. АТА... АТ- 
Bm2 е ПИЯ АТАА. ААТА, ТТ, А. ACT, А, АС ТТ. АААТАА, АТТТ. .. ААА . ATTTAA, А. АТА... -AT- 
Bm3 TEN Kers АТАА. ААТА. TT. А. AACT.A..AC TT. АААТАА. ATTT... AAA .ATTTAA. A. АТА... -AT- 
Bm4 o roam (E АТАА. ААТА, ТТ, А. ACT, A, АС ТТ. AAATAA, ATTT... AMA. ATTTAA, А. ATA.. .-AT- 
Bm5 E E ATAA. ААТА. ТТ. А. AACT. A . AC. TT. АЛАТАА. ATTT... AAA .ATTTAA. А. ATA.. .-AT- 
Bman1 TTE Ress ATAA. ААТА. ТТ. А. ААСТ. А. .АС ТТ. АААТАА. ATTT... AAA . ATTTAA, А. АТА... -AT- 
Bman2 Mendes Nis АТАА. ААТА. TT. А. АА. T. A. . AC TT. АААТАА. ATTT. А. AAA . АТТТАА. А. АТА... -AT- 
Bman3 Acus Кыш АТАА. ААТА. ТТ. А. ААСТ. А. .АС ТТ. АААТАА. ATTT... ААА . ATTTAA, А. ATA...-AT. 
П 1 900000 125:ВХАПППППППППППП 
Fig. 1 Variable sites of the alignments of the 125 rRNA gene sequences of nine species of silk- producing insects 
Dngununugn4en»rsnaunupnmmuunugaunpnsenagapnaunupnangapnagagugaug ."mnaupagnuBg -"Dü 


O 0 0 Eighty variable sites and 50 parsimony informative sites were detected from the alignments with a total of 440 bp in length. The number on the 
sequence refers to the variable sites. .[] Identical[] — 0 Indel. 
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Table З Numbers of transition/transversiod] above diagonalļ[] and genetic distance$] below diagonal|] for 12S rRNA gene sequences 





1 2 3 4 5 6 7 8 


9 10 11 12 13 14 15 16 17 18 











Ap1/2 [010 0/1 1/5 4/4 0/5 5/0 4% 6/11 10/23 9/33 10/33 9/33 8/34 9/33 9/34 9/34 8/34 34/60 
Арз [020 0.002 1/2 4/3 0/4 4/9 4/7 611 10/23 9/34 10/34 9/34 8/35 9/34 9/33 9/33 8/35 33/60 
Apr 030 0.007 0.005 5/4 2/6 4/12 4/0 7/2 11/25 11/36 9/37 10/36 10/36 10/36 10/35 10/35 9/37 32/62 
Ar [] 4[] 0.020 0.017 0.022 5/5 5/11 7/9 8/9 12/21 14/34 13/34 13/34 13/34 13/34 13/33 13/33 12/35 30/62 
Ayl 050 0.012 0.010 0.015 0.022 4/10 4/12 6/8 10/21 9/34 10/34 9/34 8/35 9/34 9/33 9/33 8/35 34/58 
Am  []6[] 0.033 0.030 0.035 0.038 0.025 2/7 4/9 7/21 8/31 8/32 8/31 8/31 8/31 8/30 8/30 8/30 31/55 
Ay2 070 0.025 0.022 0.027 0.035 0.027 0.022 4/13 8/26 7/32 8/32 7/32 7/32 7/32 7/32 7/31 6/31 30/60 
Aa 1080 0.040 0.038 0.043 0.040 0.027 0.027 0.035 8/26 8/36 9/36 8/36 8/36 8/36 8/35 8/35 7/37 32/58 
Scr [190 0.083 0.080 0.086 0.083 0.075 0.067 0.080 0.083 11/36 11/37 11/36 11/36 11/36 10/37 11/35 11/37 31/65 
Bmanl [] 100 0.108 0.111 0.117 0.119 0.105 0.097 0.097 0.111 0.114 1/0 оо оо 00 ол 0/1 1/1 31/58 
Bmi [0110 0.111 0.114 0.119 0.122 0.108 0.100 0.100 0.114 0.116 0.002 1/0 1/0 1/0 1/0 1/0 2/1 32/58 
Bm4 10 120 0.108 0.111 0.117 0.119 0.105 0.097 0.097 0.111 0.114 0.000 0.002 000 00 ол 0/1 1/1 30/59 
Bm3 0 130 0.108 0.111 0.117 0.119 0.105 0.097 0.097 0.111 0.114 0.000 0.002 0.000 0/0 ол 0/1 1/1 31/58 
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